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Abstract 

Various QCD studies based on jet observables in hadronic final states in deep-inelastic 
scattering and hadron-hadron collisions are presented. Measured quantities are event 
shape variables, jet rates and jet cross sections. QCD analyses have been performed to 
determine the value of the strong coupling constant a s (Mz), the gluon density in the 
proton and to test power corrections. 
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1 Introduction 

Hadronic final states in high energy particle collisions 
offer a rich field for testing quantum chromodynamics 
(QCD). Especially suited for such studies are infrared 
and collinear safe "jet observables" for which QCD pre- 
dictions can be calculated perturbatively. In this contri- 
butiorjLjneasuremenla are presented of event shape vari- 
ables Era, jet rates era and jet cross sections El in deep- 
inelastic scattering (DIS) by the HI and the ZEUS collab- 
oration aa well as jet cross sections in proton-antiproton 
colhsionsH by the CDF collaboration. These observables 
are directly (i.e. already at leading order) sensitive to 
the value of the strong coupling constant a s . Further- 
more the jet cross sections are directly sensitive to the 
gluon density in the proton. For all these observables, 
the perturbative cross sections have been calculated to 
next-to-leading order (NLO) in a s . Based on these data, 
QCD analyses have been performed to extract the value 
of the strong coupling constant a s (Mz) and to determine 
the gluon density in the proton. The event shape data 
have been analyzed for tests of "power corrections" . 



2 Event Shape Variables in DIS 

Event shape variables in deep-inelastic scattering events 
have been investigated over a large range of momentum 
transfers 4 < Q < 100 GeV n based on data taken at 
HERA by the HI experiment in 1994-1996 and by the 
ZEUS experiment in 1995-1997. These event shapes 
are defined by linear sums of momenta of all hadronic fi- 
nal state particles in the current hemisphere of the Breit 
frame and are in different ways sensitive to the collima- 
tion of the energy flow along the event shape axis: 
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Figure 1: The Q dependence of the mean values of 1 — Tz measured 
in deep-inelastic scattering in the current hemisphere of the Breit 
frame. The ZEUS data are measured in different bins of XBj while 
the HI results are integrated over all xej . 
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While in the analysis by the HI collaboration only the 
dependence on the four-momentum transfer Q is mea- 
sured (integrating over the Bjorken scaling variable XBj), 
the ZEUS collaboration has also studied the dependence 
on xbj. The mean values of r z are shown in Fig. | as a 
function of Q. The Q dependence of the mean values of 
the other event shapes is shown in Fig. S. It can be seen 
that at larger values of Q and XBj the final state in the 
current hemisphere of the Breit frame is more collimated 
(relative to its total energy). 
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Figure 2: Event shape variables defined in the current hemisphere 
of the Breit frame in deep-inelastic scattering. Displayed is the 
Q dependence of the mean values of the different event shapes. 
The data are compared to the results of a QCD fit involving NLO 
predictions combined with power corrections. Also shown is the 
NLO prediction. 



The QCD prediction for the mean value of any event 
shape (F) can be written as a sum of perturbative and 
non-perturbative contributions: 



(F) = (F) pcn + (F) 



(1) 



with 



(F) 



pert 



ci a s (Q) + c 2 a 2 s (Q) , 



(F) pow =a F ^| (aoQn) - a s (Q) ~ b a 2 s (Q)) . 

The perturbative part can be calculated to next-to- 
leading order. Apart from the dependence on the parton 
density functions (included in the coefficients Cj) the only 
free parameter is a s . The non-perturbative contributions 
are modeled by ap/Q, power corrections, where the 
have been calculated 0. This contribution also depends 
on the value of do(fii), the effective coupling below the 
infrared matching scale fii. 

Fits have been performed to the mean values of 
the event shapes with the free parameters a s (Mz) and 
oto(fii) according to eq. (1). The fit results in Fig. ||show 
smaller non-perturbative contributions towards higher Q. 
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Figure 3: The results for ao and a s as determined in a fit of the 
theoretical predictions based on NLO pQCD and power corrections 
to the Q dependence of the event shape mean values. 



The resulting values of the strong coupling constant 
a s (Mz) and ao (displayed in Fig. ||) are of the same 
size and demonstrate approximately universal behavior 
for all event shapes. Remaining discrepancies may be 
traced back to large theoretical uncertainties: large NLO 
corrections and a large renormalization scale dependence 
indicate possible large contributions from higher, uncal- 
culated orders. 



3 Multijet Rates in DIS 

To determine a s , the production rate of multijet events 
in DIS has been studied based on data taken in 1994 
and 1995 by the HI collaboration. Two analyses are pre- 
sented where jets are defined by the modified JADE jet 
algorithm El (mJADE), applied in the HERA laboratory 
frame. In the mJADE jet definition the final state parti- 
cles (and a pseudoparticlc to account for the non-visible 
proton remnant) are clustered in the order of increasing 
invariant di-particle masses. For this jet definition the jet 
resolution parameter is given by y^ = M 2 - / W 2 , 
where W is the total invariant mass of the hadronic fi- 
nal state and My is the smallest invariant mass between 
pairs of (pseudo-) particles. 

Two observables have been measured: the 7/2 spec- 
trumu and the (2+l)-jet ratecl, the latter as a function 
of Q 2 . The variable yi is defined as the value of the jet 
resolution parameter y^ at which exactly (2+1) jets are 
resolved (the "+1" denotes the proton remnant). The 
?/2 spectrum has been measured for momentum transfers 
200 < Q 2 < 10000 GeV 2 (Fig. @). The (2+l)-jet rate 
R2+1 is here defined as the fraction of DIS events that 
have (2+1) jets at a fixed value of the jet resolution pa- 
rameter j/ cut = 0.02 (i.e. the rate of events which have 
(2+1) jets after the clustering has been stopped when all 
Mfj/W 2 > 0.02). The results have been obtained in four 
Q 2 bins at 40 < Q 2 < 4000 GeV 2 (Fig. |top). 

Both measurements have been used to extract 
a s (Mz)- For the 2/2 spectrum, a fit has been performed to 
determine a s at a scale of fj, 2 = (Q 2 ) = 620 GeV 2 . Using 
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Figure 4: The y2 spectrum measured in deep-inelastic scattering 
with the mJADE algorithm: The data, corrected to "parton-level" , 
are compared to the NLO prediction for different values of A^ 4 ) 
using MRSH parton distributions. The best fit result is obtained for 
the value of A' 4 ) = 320 MeV corresponding to a value of a a (Mz) = 
0.118. 



the two-loop formula, this value is evolved to fi 2 . = M§ . 
For the jet rate i?2+i(Q 2 ), values of ce s have been deter- 
mined in each Q 2 -bin at a scale [i 2 , = Q 2 (Fig. [| bottom). 
A combined fit of the two-loop formula to these four val- 
ues leads to the result for ct s (Mz). The results from both 
analyses are 



y 2 : a s (M z ) 
R 2+1 : a s (M z ) 



0.118 ± 0.002(stat) iauostsyst) 



+0.007 /,v \ 
-0.006 ( theor ) 



0.117 ± 0.003(stat) t°;^(syst&theor) 
+0.006(rec. scheme) . 



The extracted a s (Mz) values are in good agreement with 
each other and with the current world average value. For 
both results the largest uncertainties come from a vari- 
ation of the parton densities used, the renormalization 
scale dependence of the NLO prediction and from a large 
model dependence of the size of the non-perturbative 
hadronization corrections. 

4 Inclusive Jet Cross Section in pp Collisions 

Inclusive jet production in pp collisions has been studied 
by the CDF collaboration based on data taken at the 
Tevatron during Run IB. The inclusive jet cross section 
has been measured □ over a wide range of transverse jet 
energies 40 < Et < 450 GeV using a cone jet algorithm 
with a cone size of R — 0.7 in the pseudo-rapidity range 
between 0.1 and 0.7. 
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Figure 5: The rate of (2+l)-jet events in in deep-inelastic scatter- 
ing as defined by the mJADE algorithm. The upper plot shows 
the data, corrected to "parton-level" together with the NLO pre- 
dictions for different values of A( 4 ) for MRSH parton distributions. 
The lower plot shows the values of the strong coupling constant 
ot s extracted in each Q 2 bin. The inner error bars show the sta- 
tistical error; the total error is given by the full error bars. The 
result of the fit of the renormalization group equation to the single 
a s (Q' 2 ) values (considering only the statistical errors) is shown as 
the dashed line, the errors as dotted lines. Taking also into ac- 
count the full systematic uncertainties results in a fit range which 
is displayed by the full lines. The central fit result corresponds to 
a value of a s (M z ) = 0.117. 
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Figure 6: The values of the strong coupling constant o: s (Et /2) and 
&s{Mz) as determined from the inclusive jet cross section in pp- 
collisions. The single values have been evolved to fi = Mz using 
the two-loop formula. The error weighted average over the Ex 
range 40 - 250 GeV gives the final result of a s (M z ) = 0.115. 
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Figure 7: The values of the strong coupling constant a a (Mz) as 
derived from a QCD fit to the inclusive jet cross section in pp- 
collisions. Several fits have been performed for different parame- 
terizations of parton density functions. Shown is the dependence 
of the central fit result of ct s (Mz) on the value of ct B (Mz) used in 
the fit of the parton distributions. 



Each of the 36 individual bins in Et allows to extract 
a value of ot s at a scale chosen to be fi = Et/2u. The 
single a s (-Er/2) values have been evolved to a s (Mz) us- 
ing the two-loop equation (Fig. |^). By forming an error 
weighted average over the Et range 40 — 250 GeV, an 
average value of a s (Mz) has been obtained (the exclu- 
sion of the points Et > 250 GeV avoids possible bias by 
the population of events at high E T ). For the CTEQ4M 
parton distributions, the result is 

a s (M z ) = 0.1152 ± 0.0001 (stat.)io;oo93(exp.syst.) 
with a theoretical uncertainty_of ±0.005 from the renor- 
malization scale dependence!!!! 

The result also depends on the parton densities used 
in the NLO calculation. This dependence is demon- 
strated in Fig. |?] where the resulting a s (Mz) value is 
shown as a function of the a s (Mz) used in the global fit 
of the parton densities. A clear correlation is seen which 
reflects the dependence of the inclusive jet cross section 
on the gluon density g{x) in the proton. The size of the 
gluon density in different global fits is anti-correlated to 
the assumed value of a s , since most cross sections only 
depend on the product a s ■ g(x). 

Due to the dependence of the extracted a s (Mz) on 
the parton densities (and their a s value) this analysis 
can not be seen as an independent measurement. How- 
ever, the fact that for standard parton distributions (as 
CTEQ4M) the extracted a s is close to the a s used in 
the parton densities is an important demonstration of 
consistency. 



5 Dijet Cross Sections in DIS 

To determine the gluon density in the proton, dijet pro- 
duction in deep-inelastic scattering has been studied hy 
the HI collaboration based on data taken in 1994-1997B. 
The dijet cross section has been measured for different 
k±-type jet algorithms in the Breit frame. QCD mod- 
els predict that hadronization corrections are smallest 
for the longitudinally boost-invariant k± jet algorithm^] 
(about 7% with negligible model dependence). 

For this algorithm the inclusive dijet cross section 
has been measured as a function of several observables: 
Q 2 , xb}, the average transverse jet energy of the two jets 
in the Breit frame Et, the invariant dijet mass Mjj and 
the variable £ = £Bj(l + M 2 j/Q 2 ) which is (at leading 
order) equal to the momentum fraction of the proton 
carried by the parton that enters the hard process. Dijet 
events have been selected where at least two jets have 
transverse energies of Et > 5 GeV in the Breit frame. 
In addition, the sum of the transverse jet energies of the 
two highest Et jets in the event has been required to be 
E T i + E T 2 > 17 GeV. 

The dijet cross section is shown in different bins of 
Q 2 as a function of Et (Fig. ||) and the variable £ (Fig. [|) 
for 10 < Q 2 < 5000 GeV 2 . The data are compared to the 
NLO prediction for different renormalization scales. Rea- 
sonable agreement is observed, except at smaller values 
of £ where the NLO prediction is slightly too low. 
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Figure 8: The double differential dijet cross sections d 2 cr/d<2 2 d&r 
in deep-inelastic scattering. The data are compared to the NLO 
prediction for two different choices of the renormalization scale for 
CTEQ4M parton distributions. 



The determination of the gluon density in a QCD fit to 
these dijet cross sections also requires knowledge of the 
quark densities and of a s . In the present analysis a s (Mz ) 
has been taken to be the world average E3 within its un- 
certainty (a s (Mz) = 0.119 ± 0.005) and subsequently 
evolved to the relevant scales using the two-loop equa- 
tion. To have a consistent treatment of the quark den- 
sities in the QCD fit, HI .data on the inclusive neutral 
current DIS cross section B at 200 < Q 2 < 650 GeV 2 
have been included in the fit. These data give strong 
constraints on the quark densities. 

The NLO corrections to the dijet cross sections at 
Q 2 < 100 GeV 2 are sizeable (50 - 100 %) and the renor- 
malization scale dependence is large (up to 20%). To 
avoid these theoretical uncertainties, only dijet data at 
Q 2 > 200 GeV 2 have been considered in the fit. In this 
range, approx. 55 % of the dijet cross section is gluon in- 
duced. To maximize the sensitivity to the x-dependence 
of the gluon density, the fit has been performed on 
the double differential cross sections d 2 a/dQ 2 d^ and 
d 2 cr/d<5 2 dxBj. The gluon and quark densities have been 
fitted at a factorization scale of [i 2 — 200 GeV 2 which 
is of the size of the hard scales for both the dijet and 
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Figure 9: The double differential dijet cross sections d 2 o-/dQ 2 d{; 
in deep-inelastic scattering. The data are compared to the NLO 
prediction for two different choices of the renormalization scale for 
CTEQ4M parton distributions. 



the inclusive cross section (dijet cross section: fi 2 ~ E T , 
inclusive cross section: /z 2 ~ Q 2 )- The ^-dependence of 
the gluon and quark densities has been parameterized by 
the usual 3-, 4- or 5-parameter functionsB. 

The gluon density obtained from the QCD fit with 
its error band is shown in Fig. ^o| in the range 0.01 < x < 
0.1. This result extends the x-range from previous gluon 
determinations from HI structure function data to larger 
x-values. The uncertainty is dominated by the error of 
a s , the renormalization scale dependence of the NLO cal- 
culation and the uncertainty in the absolute calorimetric 
energy scale. While the result is slightly higher than the 
results of different global analyses (although compatible 
within the error) it is in good agreement with the gluon 
density from a QCD analysis of HI structure function 
data. 
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Figure 10: The gluon density in the proton in the A/5-scheme 
determined in a NLO QCD fit to dijet cross sections for a value of 
a a (Mz) = 0.119±0.005. The error band is compared to the gluon 
densities obtained in global fits and to the result of a QCD fit to 
HI structure function data (only the central value without error 
band) . 
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Summary 

Measurements of jet observables in hadronic final states 
in deep-inelastic scattering and in hadron-hadron colli- 
sions have been presented. Based on perturbative QCD 
calculations in next-to-leading order, fits to the data have 
been performed to extract the value of the strong cou- 
pling constant a s (Mz) and the gluon density in the pro- 
ton. The results are consistent with each other and in 
agreement with world average values. 

The precision of these results is limited by uncertain- 
ties in the perturbative predictions (renormalization scale 
dependence) and by the correlation between the gluon 
density and a s . Further progress will require combined 
analyses of different data sets for consistent, simultane- 
ous determinations of a s and the gluon density. 

The "power correction" approach gives a fair de- 
scription of hadronization corrections to DIS event shape 
data. Further developments in this field will improve the 
understanding of non-perturbative physics and may al- 
low to describe non-perturbative contributions to differ- 
ent processes in a universal way. 
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